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A simple kinematic model based on the superposition of p + p collisions, relativistic geometry and 
hadronic rescattering is used to predict the elliptic flow observable in -y/ijvjv =2.76 TeV Pb + Pb 
collisions. A short proper time for hadronization is assumed. The predictions are compared with 
recent experimental measurements of elliptic flow which have been made for this colliding system 
and energy. It is found that the model predictions do a reasonable job in describing the experimental 
results, suggesting that the parton phase in these collisions may be short-lived. 



The CERN Large Hadron Collider has recently begun 
delivering Pb + Pb collisions at ^/snn — 2.76 TeV to 
experiments. These are the highest energy heavy-ion 
collisions ever to be produced in the laboratory. The 
LHC ALICE experiment [l| has already posted two ex- 
perimental papers based on data from these collisions, 
one in which the charged particle (hadron) multiplicity 
density is measured [2| and another in which the elliptic 
flow of charged hadrons is measured^. These are both 
important and basic observables to measure in heavy-ion 
collisions since the charged particle multiplicity density 
is related to the initial particle density and the elliptic 
flow is sensitive to the initial dynamics of the particles, 
both of which could be messengers of possible exotic phe- 
nomena taking place in these collisions [lj]. 

In order to better understand the underlying physics 
driving these observables at the LHC, a simple kinematic 
model has been constructed [U, HI with the goal of com- 
paring predictions of this model with the experimentally 
measured observables. The basis of the model is that 
the initial state of the heavy-ion collision is determined 
by the superposition of proton-proton collisions followed 
by the mutual scattering of the hadrons produced in the 
collision. Besides its simplicity, the advantages of this 
model are 1) the model has been shown to describe the 
overall trends of hadronic observables in lower energy 
Au + Au collisions at ^/s^n = 0.20 TeV from the Rela- 
tivistic Heavy Ion Collider (RHIC)[iJ, and 2) the model 
is easily scalable to LHC energies. These will be "limiting 
case scenario" predictions in the sense that only hadrons 
are used as the degrees of freedom in this model even 
at the early stages of the collision where parton (quark 
and gluon) degrees of freedom are thought to be more 
appropriate, i.e. a short proper time for hadronization 
is assumed. In spite of this assumption, it is interesting 
to note that at RHIC energies, this model has even been 
able to satisfactorily reproduce the quark number scal- 
ing seen in experiment for elliptic flow measurements for 
Au + Au collisions^]. Quark number scaling seen in ex- 
perimental observables at RHIC is often cited as evidence 
that Quark Matter is produced in RHIC collisions @ 

A prediction of the multiplicity density from this model 



has already been compared with the measured multiplic- 
ity density in central Pb + Pb collisions at ^/snn = 2.76 
TeV from ALICE in Reference 0. The outcome is that 
the model under-predicts the multiplicity density observ- 
able, dN/drj, at mid-rapidity by ~ 20% with respect to 
the lower bound of the measurement. Considering that 
the model was used to extrapolate from RHIC to LHC 
energy which represents a factor of 14 increase in y'sjvw, 
it is somewhat encouraging that the model prediction is 
even this close to the measurement, and could be consid- 
ered a validation to some degree of the simple approach 
taken in the model. 

In this paper a further and more stringent test of the 
model is presented by comparing model predictions for 
elliptic flow with the recent ALICE measurements of this 
observable [I]- The model predictions are based on the 
same calculations which were made in the prediction of 
the multiplicity density mentioned above in order to get 
a consistent picture of the ability of this simple model to 
describe the experimental results. A brief description of 
the model is presented below followed by the comparison 
of the model predictions for elliptic flow with the ALICE 
experiment for Pb + Pb collisions at y/spjN = 2.76 TeV. 

The model calculations are carried out in five main 
steps: 1) generate hadrons in p + p collisions from the 
event-generator PYTHIA, 2) superpose p + p collisions 
in the geometry of the colliding nuclei, 3) employ a sim- 
ple space-time geometry picture for the hadronization of 
the PYTHIA-generated hadrons, 4) calculate the effects 
of rescattering among the hadrons, and 5) calculate the 
hadronic observables. These steps will now be discussed 
in more detail. 

Thep + p collisions were modeled with the PYTHIA 
code [7|, version 6.409. The internal parton distribution 
functions "CTEQ 5L" (leading order) were used in these 
calculations. Events were generated in "minimum bias" 
mode, i.e. setting the low-pr cutoff for parton-parton col- 
lisions to zero (or in terms of the actual PYTHIA param- 
eter, ckin(3) = 0) and excluding elastic and diffractive 
collisions (PYTHIA parameter msel = 1). Runs were 
made at s/s = 2.76 TeV to simulate LHC collisions. In- 
formation saved from a PYTHIA run for use in the next 
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step of the procedure were the momenta and identities 
of the "direct" (i.e. redundancies removed) hadrons (all 
charge states) n, K, p, n, A, A, p, u, 77, 77', </>, and K* . 
These particles were chosen since they are the most com- 
mon hadrons produced and thus should have the greatest 
effect on the hadronic observables in these calculations. 

A main assumption of the model is that an adequate 
job of describing the heavy-ion collision can be obtained 
by superposing PYTHIA-generated p + p collisions cal- 
culated at the beam 1/i within the collision geometry of 
the colliding nuclei. Specifically, for a collision of im- 
pact parameter b, if f(b) is the fraction of the overlap 
volume of the participating parts of the nuclei such that 
f(b = 0) = 1 and f(b = 2R) = 0, where R = 1.2A 1 / 3 
and A is the mass number of the nuclei, then the number 
of p + p collisions to be superposed will be f(b)A. The 
positions of the superposed p + p pairs are randomly dis- 
tributed in the overlap volume and then projected onto 
the x — y plane which is transverse to the beam axis de- 
fined in the z-direction. The coordinates for a particular 
p + p pair are defined as x pp , y pp , and z pp = 0. The 
positions of the hadrons produced in one of these p + p 
collisions are defined with respect to the position so ob- 
tained of the superposed p + p collision (see below) . 

As was done in similar calculations for lower-energy 
RHIC collisions to give better agreement with experi- 
mental dn/dr] distributions a lower multiplicity cut 
was applied to the p+p collisions used in the present cal- 
culations which rejected the lowest 20% of the events. 
The spirit of this cut is to partially compensate for the 
fact that there is no re-interaction of primary nucleons 
from the projectile-target system in this model. 

The space-time geometry picture for hadronization 
from a superposed p+p collision located at (x pp , y pp ) con- 
sists of the emission of a PYTHIA particle from a thin 
uniform disk of radius 1 fm in the x — y plane followed by 
its hadronization which occurs in the proper time of the 
particle, r. The space-time coordinates at hadroniza- 
tion in the lab frame (xh, Hh, %hi th) f° r a particle with 
momentum coordinates (p x ,p v ,p z ), energy E, rest mass 
mo, and transverse disk coordinates (xo,yo), which are 
chosen randomly on the disk, can then be written as 
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source) has been applied to e + — e" collisions [3 • For all 
results presented in this work, r will be set to 0.1 fm/c 
as was done in applying the present model to calculating 
predictions for RHIC Au + Au collisions |5j and Tevatron 
p + p collisions [§]. 

The hadronic rescattering calculational method used 
is similar to that employed in previous studies loL 11|. 
Rescattering is simulated with a semi-classical Monte 
Carlo calculation which assumes strong binary colli- 
sions between hadrons. Relativistic kinematics is used 
throughout. The hadrons considered in the calculation 
are the most common ones: pions, kaons, nucleons and 
lambdas (71", K, N, and A), and the p, w, ij, rf , <f>, A, and 
K* resonances. For simplicity, the calculation is isospin 
averaged (e.g. no distinction is made among a 7r + , 7r°, 
and 7r~). 

The rescattering calculation finishes with the freeze 
out and decay of all particles. Starting from the initial 
stage (t = fm/c), the positions of all particles in each 
event are allowed to evolve in time in small time steps 
(At = 0.5 fm/c) according to their initial momenta. At 
each time step each particle is checked to see a) if it has 
hadronized it > th, where th is given in Eq. (4)), b) 
if it decays, and c) if it is sufficiently close to another 
particle to scatter with it. Isospin-averaged s-wave and 
p-wave cross sections for meson scattering are obtained 
from Prakash et al. 12| and other cross sections are esti- 
mated from fits to hadron scattering data in the Review 
of Particle Physics 13]. Both elastic and inelastic colli- 
sions are included. The calculation is carried out to 400 
fm/c which allows enough time for the rescattering to fin- 
ish (as a test, calculations were also carried out for longer 
times with no changes in the results). Note that when 
this cutoff time is reached, all un-decayed resonances are 
allowed to decay with their natural lifetimes and their 
projected decay positions and times are recorded. 

The rescattering calculation is described in more detail 
elsewhere [l(], 11|- The validity of the numerical methods 
used in the rescattering code have been studied using the 
subdivision method, the results of which have verified 
that the methods used are valid 
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The simplicity of this geometric picture is now clear: 
it is just an expression of causality with the assumption 
that all particles hadronize with the same proper time, 
r. A similar hadronization picture (with an initial point 



Model runs are made to be "minimum bias" by having 
the impact parameters of collisions follow the distribution 
da/db oc b, where < b < 2R. Observables are then cal- 
culated from the model in the appropriate centrality bin 
by making multiplicity cuts as normally done in experi- 
ments, as well as kinematic cuts on rapidity and px- For 
the present study, a 3200 event minimum bias run was 
made from the model for ^snn = 2.76 TeV Pb + Pb col- 
lisions which was then used to calculate hadronic observ- 
ables. As mentioned earlier, these are the same events 
which were used in the prediction of the multiplicity den- 
sity for Reference 0). 

The elliptic flow variable, V2, is defined as 



V 2 =< cos(2</>) > 



(5) 
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FIG. 1. Comparison of model V% vs. pr plots for ^Jsnn = 2.76 TeV Pb + Pb collisions with the ALICE experiment for 10 — 20% 
and 40 — 50% centrality classes. 
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where "<>" implies a sum over particles in an event 
and a sum over events and where p x and p y are the x 
and y components of the particle momentum, and x is 
in the impact parameter direction, i.e. reaction plane 
direction, and y is in the direction perpendicular to the 
reaction plane. The V 2 variable is calculated from the 
model using Eq. ([5]) and taking the reaction plane to be 
the model x — z plane. As seen, if < p x >^< p y >, then 
V 2 ~ 0, and for < p x »< p y >, then V 2 — 1. 

Figure [T] shows the comparison of model V2 vs. pt plots 
for y/sNN = 2.76 TeV Pb + Pb collisions with measure- 
ments from the ALICE experiment^ for 10 — 20% and 
40 — 50% centrality classes. Note that all model predic- 
tions in this paper are compared with the experimental 
value, which is V 2 extracted using the 4-particle cu- 
mulant method since this observable minimizes non- 
flow correlations in the elliptic flow which are not present 
in the model, thus making it more comparable. The 
model predictions are seen to follow the experimental 
points fairly closely within the uncertainties shown for 
both centrality classes. The uncertainties on the model 
prediction are statistical. The model accurately describes 
the increasing V 2 with increasing px for px < 2 GeV/c 
and then the "flattening" of the dependence of V 2 on fur- 
ther increase of pr for Pt > 2 GeV/c, although somewhat 
under-predicting the measurement for px > 2 GeV/c for 
the 10 — 20% centrality class. From the model picture all 
of these dependences, as well as the non- vanishing values 
of V 2 , are a result of the hadronic rescattering. If rescat- 



tering is turned off in the model, or equivalently r ^> 0.1 
fm/c, V 2 -> for all cases 0, [1]. 

Figure [2] compares V 2 integrated over 0.2 < pr < 
5.0 GeV/c vs. centrality class with the ALICE 
experiment 0. The model is seen to qualitatively de- 
scribe the trend of the measurements of increasing V 2 
with increasing centrality class and then reaching a max- 
imum, but to under-predict the measurements on average 
by ~ 14%. This under-prediction could be a reflection of 
the under-prediction of the model for the particle mul- 
tiplicity, mentioned earlier (which could be corrected in 
the model by including some degree of re-interaction of 
the primary nucleons from the projectile-target system). 

From the comparisons shown above in Figures [1] and 
[5] between the model and experiment, it is seen that the 
model does a reasonably good job in describing the main 
features of the experimental elliptic flow results for LHC 
Pb + Pb collisions. The model even shows the "high-py 
flattening" feature of V 2 which pure hydrodynamic mod- 
els that are based on the formation of Quark Matter in 
these collisions and thus use parton degrees of freedom in 
the early stage of the collision are not able to reproduce 
[lij . While it is interesting and perhaps a little disturb- 
ing that so simple a model is able to do as well as it 
does in describing the experimental results, the question 
of "what physics have we learned from this?" could be 
posed. The main motivation for this model has been to 
study to what extent hadronic observables in relativistic 
heavy-ion collisions are influenced by hadronic rescatter- 
ing. It is expected that some hadronic rescattering will 
take place in these collisions after some transition from 
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FIG. 2. Comparison of model integrated Vn vs. centrality 
class for yi^T = 2.76 TeV Pb+Pb collisions with the ALICE 
experiment. 



hadronization occurs. Further model studies and com- 
parisons with experiment for other observables may help 
clarify which of these interpretations is the more valid. 

In conclusion, a kinematic model based on the super- 
position of PYTHIA-generated p + p collisions, relativis- 
tic geometry and hadronic rescattering has been used in 
the present work to predict the elliptic flow observable in 
y/SNN = 2.76 TeV Pb + Pb collisions. A short proper 
time for hadronization of t = 0.1 fm/c has been as- 
sumed as in previous studies with this model which have 
shown qualitative agreement with experiments. It has 
been shown that the simple hadronic rescattering model 
accurately describes the features of the ALICE elliptic 
flow measurements, suggesting that hadronic rescatter- 
ing plays an important role in determining the properties 
of the elliptic flow observable in these collisions. These 
results also suggest that the parton phase in these colli- 
sions may be short-lived. 
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a possible Quark Matter state to hadronic matter. As 
mentioned earlier, the "limiting case scenario" has been 
studied in the present model where only hadronic degrees 
of freedom are considered. An initial-state model is used 
to set up the initial kinematics, in the present case a su- 
perposition of PYTHIA p + p collisions due to its conve- 
nience but which alternatively could have been a thermal 
or parton model as have been used in a past successful el- 
liptic flow model study for RHIC collisions 0, but the 
hadronic rescattering is the "active ingredient" driving 
the dynamics of the hadronic observables in the model. 
In order to get enough rescattering to agree with experi- 
ment, the rescattering must begin almost immediately in 
the collision. This suggests three possible interpretations 
as to the meaning of the model at the early stage of the 
collision where its assumptions are the most contentious: 
1) the hadronization time is very short, i.e. r = 0.1 fm/c, 
and thus hadronic degrees of freedom dominate from the 
beginning, 2) the early stage is partonic and more ex- 
tended in time and the early-stage hadronic scattering 
in the model "mocks up" some qualitative features of 
the early-stage parton scattering (including "viscosity" 
in the language of hydrodynamic models), or 3) a com- 
bination of 1) and 2) wherein a mixed phase of partons 
and hadron-like objects initially coexist until complete 
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